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Abstract 


Experimental models for estimating erosion have been developed for a specific area and their 
calibration is necessary for use in conditions other than their own location. Examining the accuracy of 
experimental models for estimating erosion can lead to better estimates of sediment load and thus 
better design of soil and water protection operations. Therefore, it is necessary to identify high-risk 
areas of erosion to control and reduce erosion and sediment production. The study aimed to investigate 
the accuracy and capability of ICONA, support vector machine, Chaid and random forest models in 
estimating erosion. First, digital layers of variables affecting erosion including slope, geological 
formation, land use, soil, height, slope direction, surface curvature, waterway network density, 
distance from the waterway, fault density, distance from the fault, and topographic moisture index 
(Twi) were prepared. To compare different models, statistical indices of correlation coefficient (R) 
and absolute magnitude of error (MAE) were used. The results showed that among the mentioned 
models, the support vector machine model, ICONA and random forest with M7, M9 and M12 pattern 
had the highest accuracy with correlation coefficient of 0.899, 0.845, and 0.921 and the lowest mean 
absolute value. It has error MAE = 0.711, MAE = 0.721, and MAE = 0.628. According to the study 
of effective factors in soil erosion model, it is concluded that the parameters of the slope, geological 
formation, land use, soil, distance from the waterway, and topographic moisture index (Twi) are more 
sensitive to erosion and the factors affecting erosion in these areas are more active. Most of the study 
area is part of a very high to high erosion class that these classes are mainly located in the center of 
the area. Most areas at high to severe erosion risk are located in the sloping topographic unit. 
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حکیده 


مدل‌های تجربی برآورد فرسایش برای یک منطقه خحاص agi‏ شده‌اند و واسنجی آن‌ها برای استفاده در 
شرایطی غیر از محل تهیه آن‌ها ضروری است. بررسی دقت مدل‌های تجربی برای تخمین فرسایش می- 
تواند منجر به برآوردهای بهتر از بار رسوبی و درنتیجه طراحی بهتر عملیات حفاظت خاک و آب گردد. 
لذا شناسایی مناطق پر خطر فرسایش جهت JES‏ و کاهش فرسایش و تولید رسوب ضروری است. هدف 
در تخمین فرسایش است. ابتدا لایه‌های رقومی متغیرهای موثر در فرسایش شامل شیب. سازند 
زمین شناسی» کاربری اراضی. خاک. ارتفاع» جهت شیب انحنا سطح» تراکم شبکه آبراهه فاصله از آبراهه» 
مقایسه مدل‌های مختلف از شاخص‌های آماری ضریب همبستگی R)‏ و قدر مطلق خطا (MAE)‏ استفاده 
شده است. بررسی نتایج نشان دادند که از Ole‏ مدل‌های نام برده مدل ماشین بردار پشتیبان» ICONA‏ 3 
جنگل تصادفی با الگوی 7 M9‏ و M12‏ بیشترین دقت را با ضریب همبستگی (AA)‏ (۸۶۵/) و 
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)0/411( و کمترین ميانگين قدر مطلق خطا (۰/۷۱۱= MAE <۰/۷۲۱( (MAE‏ و (MAE =:/1vA)‏ 
دارد. با توجه به مطالعه عوامل مؤثر در مدل فرسایش خاک این نتیجه به دست می‌آید که پارامترهای شیب» 
سازند زمین‌شناسی» کاربری اراضی» SLE‏ فاصله از آبراهه و شاخص رطوبت توپوگرافی (Twi)‏ 
حساسیت بیشتری نسبت به فرسایش داشته و عوامل موثر بر فرسایش در این مناطق فعالیت بیشتری دارند. 
بخش اعظم منطقه موردمطالعه. جزو طبقه فرسایشی خیلی زياد تا زياد قرار دارد که به‌طور عمده این 
کلاس‌ها در مرکز منطقه قرار گرفته‌اند. بیشترین مناطق تحت خطر فرسایش بالا تا شدید. در واحد 
توپوگرافی با شیب زیاد واقع شده‌اند. 


کلیدواژه‌ها: فرسایش خاک» حوضه گناباده مدل‌های تجربی» شاخص‌های آماری» ضریب همبستگی. 


Andie -۱ 


فرسایش SE‏ یکی از عوامل مهم تخریب خاک در دنا است که اثرات منفی شدیدی روی SLE‏ ها دارد. این 
از نظر وسعت و شدت فرسایش خاک در dle‏ کشورهای بسیار سیب دیده جهان قرار دارد (کرمی و همکاران 
hs aA‏ تاه تین یی اما هماع انا وراد شمه سای مان 
خاک در مناظر می‌شود (بورلی و همکاران ؛ ۲۰۲۱). در دهه‌های اخیر فرسایش خاک توسط آب با تغییر آب‌وهوا 
و کاهش تدریجی منابع طبیعی نسبت به جمعیت‌های انسانی به یک مسئله جهانی تبدیل شده است. علاوه بر این 
خاک به شیوه‌های مختلف در معرض خطرات بیشتری ازجمله فعالیت‌های نامناسب کشاورزی, قطع درختان جنگلی» 
است و از طرفی داشتن اطلاعات مکانی در ارتباط با پدیده‌ها این امکان را فراهم می‌کند که همبستگی بین پدیده‌ها و 
الگوی فضایی بین ol‏ ها شناسایی شود. درنتیجه شناخت الگوی فضایی فرسایش خاک یک WS‏ برای درک دینامیک 
و شناسایی دلایل زیستی و غیر زیستی فرسایش خاک است و تعیین الگوی پراکنش شدت فرسایش و رسوب می- 
تواند به ارائه راهکارهای مناسب برای کاهش و مقابله با فرسایش و رسوب در حوزه آبخیز کمک کند (مزبانی و 
همکاران۱۶۰۰). استفاده از مدل‌های فرسایشی, رویکردی بهینه در زمینه بررسی فرسایش, تولید رسوب و شناسایی 
مناطق بحرانی محسوب می‌شود (رضازاده و همکاران» (AYAV‏ در مدل‌های مختلف فرسایش «SL‏ شاخص‌های 
متفاوتی برای فرسایش‌پذیری توسعه oala‏ شده است. استفاده از مدل‌های فرسایشی» رویکردی بهینه در زمینه بررسی 
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پیش‌بینی برای ارزیابی هدر رفت SE‏ برنامه‌ریزی حفاظت. موجودی فرسایش خاک و برنامه‌ریزی پروژه‌ها 
مورداستفاده قرار گیرند. مدل‌ها می‌توانند به‌عنوان ابزاری برای درک فرایندهای فرسایش و تأثیر آن‌ها مورداستفاده قرار 
گیرند (موشه و همکاران" ۲۰۱۳). در این راستا مطالعات متعددی در داخل و خارج کشور انجام گرفته که به چند 
مورد از آن‌ها اشاره می‌شود. کریستو و همکاران (۲۰۱۶) بر اساس ویژگی‌های زمین فضایی مدل‌های برآورد 
فرسایش. آن‌ها را طبقه‌بندی نمودند و نتیجه گرفتند که برای دست کاری مقیاس جهت افزايش کارایی مدل تخمین 
فرسایش» دو نظریه فضایی سلسله مراتبی و فراکتال برای بهره‌برداری بیشتر در آینده پيشنهاد می‌شود. ترن تی ترن 
تی فونگ و همکاران (۲۰۱۶) برای شبیه‌سازی فرسایش خاک بالا دست حوضه رود بو در ویتنام مرکزی» بریک و 
همکاران TV‏ برای برآورد بار رسوب در حوضه کالایا در شمال کشور مراکش, دورو و همکاران*(۲۰۱۸) برای 
پیش‌بینی جریان رودخانه و بار رسوب حوضه LIST‏ در ترکیه از مدل SWAT‏ استفاده کردند. نتایج مطالعات OLS‏ 
دادند این مدل در شبیه‌سازی هیدرولوژیکی و برآورد بار رسوب می‌تواند در حوضه‌ای با ویژگی‌های مشابه 
مورداستفاده قرار گیرد. یمانی و همکاران (۱۳۹۷) در مقاله‌ای تحت عنوان کارایی SLT‏ کمی پارامترهای 
ژئومورفومتریک در تهیه نقشه حساسیت فرسایش خاک به این نتایج دست يافتند که مدل ۷150۸دارای دقت پیش - 
بینی بالاتری به روش فاکتور ترکیب است. آرشاد اشرف*(۲۰۲۰) در پژوهشی تحت عنوان مدل‌سازی حطر فرسایش 
خاک در شرایط مختلف کاربری اراضی و بارندگی در حوضه رودخانه سان به این نتایج دست یافت که عوامل 
تأثیرگذار منجر به نرخ فرسایش بالاتر تحت سناریوهای متغیر تغییر محیطی است. شدت فرسایش در حدود ۱۵ تن 
در هکتار در سال در خاک آزاد و ۱۰/۳ تن در هکتار در سال در اراضی کشاورزی پیش‌بینی شده است. کومار و 
همکاران" (۲۰۱۹) در پژوهشی تحت عنوان بررسی مقایسه‌ای مدل های مختلف فرسایش خاک و عملکرد رسوب در 
حوضه چاتیسگر هند به این نتایج دست يافتند که مدل RUSLE‏ برای فرسایش خاک و مدل MUSLE‏ برای 
عملکرد رسوب قابل‌اطمینان‌تر است. با در نظرگرفتن لزوم مطالعه فرسایش و تخمین مقادیر هدر رفت خاک توسط 
عوامل تأثیرگذار فرسایندگی و فرسایش‌پذیری» این پژوهش سعی دارد تا با استفاده از مدل‌های ICONA‏ ماشین 
بردار پشتیبان چاید و جنگل تصادفی مقادیر پارامترهای Fhe‏ در فرسایش خاک را برآورد نماید و در پایان با تلفیق 
لایه‌های رقومی عوامل موثر در نقشه پهنه‌بندی فرسایشی ag‏ شود؛ بنابراین فرضیات پژوهش حاضر ASE‏ بیشتر بر 
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۲مواد و روش‌ها 

در پژوهش حاضر برای ارزیابی فرسایش خاک حوضه از چهار مدل 100۷۸ ماشین بردار پشتیبان چاید و 
جنگل تصادفی استفاده گردید و بین نتایج این مدل‌ها در برآورد فرسایش خاک مقایسه انجام گرفت. داده‌های ارتفاع» 
شیب جهت شیب و انحنا سطح از مدل رقومی ارتفاع با قدرت SSH‏ ۲۰ متس داده‌های تراکم شبکه آبراهه» فاصله 
از آبراهه» تراکم گسل و فاصله از گسل از نقشه توپوگرافی ۸/۵۰۰۰۰ نقشه زمین شناسی ۰۱/۲٥۰۰۰‏ کاربری اراضی از 
تصاویر ماهواره‌ای لندست ETM?‏ (۲۰۱۱-۲۰۲۰) با دقت مکانی ۱۵ متر. سازندهای لیتولوژی از نقشه زمین شناسی 
۰ خاک از نقشه خاک شناسی ۱/۱۰۰۰۰۰ استخراج شدند. به منظور نمره دهی به عوامل فرسایش در این 
چهار مدل از لایه‌های اطلاعاتی شیب سازندهای زمین شناسی» کاربری اراضی» خاک ارتفاع. جهعت شیب انحنا 
سطح» تراکم شبکه آبراهه» فاصله از آبراهه» تراکم گسل» فاصله از گسل و رطوبت توپوگرافی Twi)‏ استفاده گردید. 

۱-۲- موقعیت منطقه موردمطالعه 

شهرستان گناباد در جنوب استان خراسان رضوی و در محدوده طول جغرافیایی ۵۸ درجه و عرض جغرافیایی 
۶ درجه و ۳۰ adda‏ قرار دارد. مساحت این شهرستان ۵٩۰۲‏ کیلومترمربع است. این شهرستان در اقلیم خشک و 
نیمه‌حشک حاشیه کویر قرار دارد (شکل ۱). واحدهای سنگی منطقه شامل سنگ‌های آذر آواری همانند ایگنمبریت» 
آگلومره توف اسیدی» سنگ‌های آتشفشانی شامل ریولیت» داسیت و تراکیت (سنوزوئیک) و رخساره‌های شيل و 
ماسه سنگ شمشک (ژوراسیک میانی) است. از نظر ساختاری ناحیه مورد نظر با توجه به تقسیم‌بندی پهنه‌های 
رسوبی- ساختاری عمده ایران در شمال شرقی خرده قاره ایران مرکزی قرار دارد. سازند شمشک (ژوراسیک میانی) 
بیشترین گسترش را در سطح منطقه دارد. ساختمان‌های زمین‌شناسی متنوع مانند چین‌های برگشته چین‌های 


خوابیده گسل‌های راست لغز و راندگی در ناحبه یافت می‌شود. 
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شکل ۱- موقعیت جغرافیایی منطقه 

۲-۲- روش انجام پژوهش 

ICONA Ju ۱-۲-۲ 

مدل مذکور توسط انجمن علمی حفاظت از طبیعت اسپانیا ارائه شده است. این مدل در کشورهای اروپایی و 
بسباری از مناطق مدیترانه‌ای قابل استفاده بوده و در Ole‏ بسباری از راه‌های موثر بر پیش‌بینی فرسایش با استفاده از 
پذیرفته شده است (آی کونا: ۱ این مدل شامل هفت مرحله است. مرحله یک: لابه شیب. مرحله دو: لایه 
گیاهی» مرحله شش: لایه حفاظت خاک و لابه هفتم: لایه فرسایش آبی است. در ابتدا نقشه شیب و نقشه 
زمین‌شناسی تهیه می گردد و در ادامه از همپوشانی لایه شیب و زمین‌شناسی لایه فرسایش‌پذیری خاک تهیه می‌شود. 
در مراحل بعدی نقشه کاربری و پوشش گیاهی تهیه می‌شود و با همپوشانی لابه کاربری اراضی و پوشش ALS‏ 
نقشه حفاظت خاک به دست می‌آید. سپس AY‏ فرسایش پذیری خاک و a‏ حفاظت خاک همپوشانی شده و نقشه 


1 ICONA 
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۲-۲-۲- مدل ماشین بردار پشتیبان 

ماشین بردار پشتیبان یک سیستم یادگیری کارآمد بر مبنای تثوری بهینه‌سازی مقید است که از اصول استقرای 
کمینه‌سازی lee‏ ساختاری استفاده کرده و منجر به یک جواب بهینه IS‏ می‌شود. در یک مدل رگرسیونی SVM‏ 
لازم است وابستگی تابع متغیر وابسته (Y)‏ به مجموعه‌ای از متغیرهای مستقل () تخمین زده شود. فرض براین است 
که مانند So‏ مسائل رگرسیونی؛ رابطه بدین متغیرهای وابسته و مستقل توسط یک تابع معين به‌علاوه یک مقدار 
اضافی نویز" مشخص می‌شود. 

y= f(x) + noise A رابطه‎ 

موضوع اصلی پیدا کردن فرم تابع ٤‏ است که بتواند Sy pod‏ صحیح موارد جدیدی را که SVM‏ تا کنون تجربه 
نکرده است پیش‌بینی AS‏ در روش ESVM‏ از توابع هسته‌ای (Ob S)‏ به‌منظور جداسازی صفحات و بهینه‌سازی 
فواصل موجود بین داده‌های مشاهداتی و تخمینی استفاده می‌شود اما در مدل SVM‏ -۷ از روش برنامه‌ریزی غير 
b>‏ برای کاهش خطای پیش‌بینی استفاده می‌شود (سیاسر و همکاران: ۱۳۹۸). در این مطالعه مدل 6-9۷۷۲ به دلیل 


کاربرد گسترده آن در مسائل رگرسیونی استفاده گردید. برای این مدل تابع خطا ب‌صورت رابطه ۲ تعریف می‌شود: 


C= YL, i * + cD is Ge + ww ۲ رابطه‎ 

تابع حطای فوق لازم است با توجه به محدودیت‌های زیر حداقل گردد: 

p(Xi)+b-yi <e +éi Y رابطه‎ 
wt 

yi-W! P(Xi)+ b Set éi :۶ رابطه‎ 

či.čix>0 :۵ رابطه‎ 


روابط 6 ثابت گنجایش» W‏ بردار ضریب» WT‏ ترانهاده بردار ضرایب. Gigi‏ ضرایب کمبود 9 ضریب ابت. N‏ 
الگوی آموزش مدل» ۵ تابع BS‏ است. 

ماشین‌های بردار پشتیبان برای حل مسائل غیرخطی. ابعاد مسئله را از GLb‏ توابع کرنل تغییر می‌دهند. انتخاب 
کرنل برای SVM‏ به حجم داده‌های آموزشی و ابعاد بردار Shy‏ بستگی دارد؛ به‌عبارت‌دیگر باید با توجه به این 
پارامترها تابع کرنلی را انتخاب نمود که توانایی آموزش برای ورودی‌های مسئله را داشته باشد. در عمل LES‏ خطی 


به کار برده می‌شود. در حدول A)‏ معادلات or‏ از کرنل‌های رایج ارائه شد: 


1 Noise 


حدول ۱- توابع کرنل رایج در ماشین‌های بردار پشتیبان (کیزی و همکاران: ۲۰۱۵) 
نوع تلع تابع ES‏ 
خحطی Xj) = x Xj‏ .106 


درنهایت. تابع تصمیم رگرسیون بردار پشتیبان غیرحطی به‌صورت رابطه 1 خواهد بود که کنترل کننده ميزان 
نوسان تابع گوسی همچنین کنترل کننده نتایج پیش‌بینی و تعمیم دهنده مدل SVM‏ است. 

f(x) <<) di — di *)k (Xi. Xj) +b A رابطه‎ 

۳-۲-۲- مدل جنگل تصادفی 

مدل جنگل تصادفی بر خلاف مدل‌های کلاسیک چون رگرسیون که تنها بر پایه یک مدل تکیه دارند با استفاده 

از صدها و هزاران درخت از اطلاعات پیشتری در داده‌ها استفاده می‌کند تا بتوان استنباط بهتری از متغیرها داشت. این 
روش یک تکنیک مدل ناپارامتری" و متعلق به خانواده روش‌های دسته‌جمعی است. این الگوریتم ازجمله 
دسته‌بندی‌هایی است که متد Bagging‏ را به کار می‌گیرد. در طی فرآیند ساخحت درخت در هر شاخه از بین تمام M‏ 
متغیر مستقل به‌صورت تصادفی ۳۵ متغیر برای تقسیم شدن انتخاب می‌شود. برای حالت رگرسیونی نسبت او با 
یک سوم است و برای کلاسه‌بندی برابر با mE VM‏ پیشنهاد شده است. روش درخت رگرسیون جنگل تصادفی 
به‌ویژه هنگامی که تعداد مشاهدات در مقایسه با تعداد پیش‌بینی کننده‌ها نسبتاً کم باشد یک روش پیش‌بینی کارآمد 
است (پهلوان و همکاران: ۲۰۱۶). محاسبات مدل جنگل تصادفی در نرم‌افزار Matlab2018‏ انجام شد. به‌منظور 
انتخاب مناسب‌ترین خط برازش یافته توسط این روش در احتمالات متفاوت. متغیر اندازه گره که نشان دهنده تعداد 
برگ‌ها در هر شاخه است با آزمون‌وخطا تعیین شد. 

CHAID الگوریتم‎ -۶-۲-۲ 

چاید مخفف مجذور مربع شناسایی کننده تعاملات خودکار است. این روش آماری بسیار موثر برای بخش‌بندی 
و رشد درخت. به‌وسیله کاس توسعه یافته است. با استفاده از مفهوم تست آماری به‌عنوان یک line‏ چاید تمام 
مقادیر صفت پیش‌بینی کننده بالقوه را ارزیابی می‌کند. این الگوریتم مقادیری را که به‌صورت آماری همگن در نظر 
گرفته شده‌اند (شباهت) با توجه به متغیر هدف و حفظ plod‏ مقادیری که ناهمگن(غیرمشابه) هستند ادغام می‌کند. 
سپس این الگوریتم بهترین پیش‌بینی کننده را از شاخه در درخت تصمیم انتخاب می‌کند به‌طوری‌که هر فرزند گره از 
یک گروه از ارزش‌های همگون نسبت به صفت انتخاب شده می‌باشند. این فرآیند به‌صورت بازگشتی تا آنجا که 
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2 Nonlinear 
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درخت به‌طور کامل رشد می‌کند. ادامه می‌یابد. آزمون آماری مورداستفاده بستگی به سطح اندازه‌گیری میدان هدف 
دارد. اگر صفت هدف پیوسته باشد. یک آزمون ۴ صورت می‌گیرد. اگر صفت هدف دسته‌ای باشد. یک آزمون 
مجذور مربع به کار برده می‌شود. چاید یک روش باینری نیست و می‌تواند بیش از دو دسته در هر سطح خاص از 
درخت ایجاد کند؛ بنابراین گرایش به ساخت درخت نسبت به روش رشد باینری گسترده‌تر است. این الگوریتم برای 
تمام انواع متغیرها پاسخگو است و دو متغیر وزنی و فراوانی را می‌پذیرد. مدل مذکور همچنین مقادیر مفقودی را با 
یک طبقه‌بندی معتبر ایجاد می‌کند (کیزی و همکاران : ۲۰۱۵). 


۳- نتایج و بحث 


۱-۳- بررسی عوامل مؤثر در وقوع فرسایش خاک 

تراکم شبکه آبراهه: هر چه قدر تراکم شبکه زهکش phy‏ شود به میزان فرسایش خاک هم اضافه می‌شود. 
کمترین میزان فرسایش خاک در پهنه‌هایی واقع شده‌اند که تراکم شبکه زهکش کمتر بوده است. تراکم زهکشی رابطه 
معکوسی با میزان نفوذپذیری دارد. تراک بالای شبکه زهکشی Shy‏ تمرکز آب به‌صورت رواناب و کاهش نفوذ 
است. از طرفی می‌توان گفت که وجود تراکم بالای شبکه زهکشی می‌تواند نشان دهنده سازندهایی با استعداد ایجاد 
جریان‌های متعدد باشد. می‌توان این‌طور استنباط نمود که علت فرسایش خاک با درجه زياد عامل تراکم شبکه 
زهکش است. بیشترین میزان خطر فرسایش خاک در تراکم زهکش با مقدار تراکم از 4۷۸۸ متر تا ۷۲۲۳ متری شبکه 
زهکش رخ داده است (شکل ۲). 

خاک: مقاومت خاک نسبت به فرسایش تا اندازه‌ای به وضعیت پستی‌وبلندی, درجه شیب و اثر تخریبی انسان‌ها 
بستگی دارد؛ اما فرسایش‌پذیری به عواملی از قبیل Cab‏ خاک پایداری دانه خاکدانه‌هاه مقاومت برشی» ظرفیت نفوذ 
و مقدار مواد آلی و شیمیایی نیز بستگی دارد (مزبانی و همکاران ۱1۰۰: (EV‏ حاک‌های آنتی سول در تشکیلات 
زمین‌شناسی جدید به‌ویژه در شیب تند که فرسایش شدید بوده یافت می‌گردد. خاک‌های آنتی سول در منطقه 
موردمطالعه در روی کوه‌ها و تپه‌ها به‌صورت خاک خیلی کم عمق يا کم عمق وجود دارد. این نوع خاک. خاک‌های 
با نفوذپذیری زیاد و هرز آب خیلی کم دارای زهکشی طبیعی مناسبی بوده و تراوش پذیری آن‌ها حتی هنگامی که 
رطوبت زیاد است و بنابراین پتانسیل تولید رواناب کمی دارند و خاک‌های با بافت سبک. شنی و دارای سنگریزه و 
قلوه‌سنگی عمیق تا خیلی عمیق را شامل می‌شوند. نفوذپذیری این خاک‌ها بین ۷/۵ تا ۱۱/۸ سانتی‌متر در ساعت متغیر 
است. خاک هیستوسول توسط بادهای خشک باعث ایجاد فرسایش می‌شود. خاک هایی هستند که روان آب متوسطی 
تولید می‌کنند. این خاک‌ها هنگامی که SUIS‏ خیس باشند قابلیت نفوذ متوسطی دارند. در این گروه خاک‌ها می‌توان 


خاک‌های با بافت نسبتاً ریز مانند خاک‌های شنی لومی خاک‌های شنی همراه با مقدار کمی رس را نام برد. نفوذپذیری 
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سال دهم بررسی یتانسیل مناطق مستعد فرسایش با مدل‌های .... ۱۱ 


آن‌ها بین ۳/۸ تا ۷/۵ سانتی‌متر در ساعت است. خاک اینسپتی سول جزو SLE‏ جوان و SLE‏ مخصوص مناطق 
خشک است که به شدت تحت تابر فرسایش قرار دارد. خاک هانی هستند که آب دوی نسبتاً زیادی تولید می‌کنند. 
قابلیت نفوذ این خاک ها هنگامی که Wal‏ خیس باشند نسبتاً کم است. اغلب دارای بافت سنگین تا خیلی سنگین می 
باشند؛ مانند خاک‌های مخلوط رس و سیلت و همچنین خاک‌های لوم رس است. در پروفیل برخی از این خاک‌ها 
یکلایه غیرقابل نفوذ وجود دارد. نفوذپذیری این خاک‌هابین ۱/۳ تا ۳/۸ سانتی‌متر در ساعت است. خاک 
ورتی سول شرایط مناسبی را جهت فرسایش ایجاد می‌کند؛ سار این نیازمند فرایند مدیریتی eke‏ است تا 
از فرسایش جلوگیری شود. خاک هایی هستند که روان آب زیادی تولید می کنند و قابلیت نفوذ آن‌ها هنگامی که 
مرطوب باشند کم است. نفوذپذیری آن‌ها کمتر از ۱/۳ سانتی‌متر در ساعت است. این گروه عمدتاً شامل خاک‌های 
رس با ظرفیت آماس‌پذیری obj‏ خاک‌هایی که عمق آب زیرزمینی در آن‌ها دائماً بالاست. خاک‌هایی که در آن‌ها لاه 
رس سخت نزدیک سطح خاک است و همچنین خاک‌هایی کم عمق بر روی مواد غیرقابل نفوذ می‌باشند. خاک 
اکسی سول خاک‌های شدیداً هوازده است که در این نوع خاک خطر فرسایش خاک به‌وسیله آب بیشتر است. خلل 
و فرج آن کم تا نسبتاً زیاد است. دارای ۳۵-6۵ درصد سنگریزه ریزودرشت و قطعات سنگ است و نفوذپذیری آن‌ها 
کمتر از ۳/۸ سانتی‌متر در ساعت است (شکل AY‏ 

شیب: این عامل بیانگر تأثیر عوامل ناهمواری و [RE‏ زمین در فرسایش خاک است. ميزان خاک ازدست‌رفته 
ناشی از عوامل ناهمواری و شکل زمین نظیر شیب و طول دامنه است. شیب زیاد باعث افزايش سرعت آب 
رودخانه‌ها و شعبات Of‏ می‌شود که این امر منجر به فرسایش‌پذیری زیاد خاک‌ها و رسوبات می‌شود. بیشترین میزان 
فرسایش خاک در شیب بیش از ۲۵ درجه رخ داده است (شکل ۲). 

زمین‌شناسی: با توجه به اهمیت ویژگی‌های سنگ‌شناسی در فرسایش و تولید رسوب بررسی حساسیت واحد- 
های سنگ و خاک به فرسایش از اهمیت بالقوه‌ای برخوردار است. سنگ‌ها رفتار متفاوتی در برابر هموازدگی و 
فرسایش از خود نشان می‌دهند. بیشترین میزان فرسایش در واحدهای زمین‌شناسی سنگ‌های آذرین شیل و ماسه 
سنگ مشاهده شده است. علت فرسایش‌پذیری زیاد خاک به دلیل سست بودن رسوبات است؛ اگرچه بیشترین واحد 
زمین‌شناسی را Ghee‏ رسی تشکیل می‌دهد. ولی به دلیل نفوذپذیری زياد فرسایش کمتر صورت گرفته است 
ARS)‏ 

کاربری اراضی: پوشش گیاهی یکی از پارامترهای مؤثر در میزان فرسایش است. با افزایش پوشش ALS‏ میزان 
فرسایش نیز کاهش خواهد یافت. چنانچه منطقه پوشیده از گیاهان انبوه فعالیت‌های کشاورزی معمول نباشد و 
سطح Ol‏ کمتر مورد تعلیف دام قرار گیرده AB‏ عامل GB‏ در میزان فرسایش خاک و تولید رسوب منفی بوده و از 
میزان فرسایش خاک و تولید رسوب عوامل دیگر می‌کاهد ولی برعکس چنانچه در بیش از نیمی از سطح یک 
منطقه, عملیات کشاورزی بدون رعایت Spel‏ حفاظت آب‌وخاک صورت گیرد با چرای سنگینی اعمال شود نقفش 
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عامل فوق در ایجاد فرسایش خاک و تولید رسوب مثبت بوده و باعث فرسایش و تولید رسوب می‌شود. اراضی 
فاقد پوشش JE‏ زیادی در فرسایش خاک خواهند داشت. چون سطح زمین فاقد هرگونه پوشش گیاهی در برابر 
عوامل فرسایش آبی و بادی است. با افزایش مساحت اراضی فاقد پوشش CALS‏ زراعت آبیء باغات و اراضی 
مسکونی فرسایش به علت کاهش نفوذپذیری افزایش پیدا می‌کند. بیشترین میزان فرسایش در مراتع فقیر به دلیل نبود 
پوشش گیاهی رخ داده است (شکل AY‏ 

ارتفاع: بیشترین میزان فرسایش در ارتفاعات بیشتر از ۲۰۰۰ متر رخداده است. علت این فرسایش‌پذیری بسیار 
زیاد تا زياد به دلیل تأثیری است که ارتفاع برافزایش بارندگی دارد. با افزايش ارتفاع بارندگی بیشتر می‌شود و بارندگی 
بیشتر به مقدار جریان Ol‏ رودخانه‌ها و شعبات آن می‌افزاید درنتیجه قدرت فرسایندگی رودخانه‌ها زیادتر می‌شود و 
خاک‌ها فرسایش می‌پابند. همچنین در ارتفاعات زیاد به دلیل عمل انقباض و انبساطی که در سنگ‌ها و خاک‌ها رخ 
می‌دهد. باعث خرد شدن و سست شدن این سنگ‌ها و خاک‌ها می‌شود و درنتیجه مواد و مصالح لازم برای حمل 
توسط رودخانه‌ها فراهم می‌آید (شکل AY‏ 

جهت شیب: شیب‌های شمالی کمتر در معرض GLU‏ خورشید بوده و درنتیجه OLALS‏ با رطوبت کافی رشد 
نموده و تبخیر و تعرق کم می‌شود ولی عمل فرسایش در شیب‌های جنوبی که آفتاب‌گیر است بیشتر دیده می‌شود 
(شکل ۲). 

انحنا: انحنا نشان دهنده میزان انحراف سطح از صاف بودن و یا به‌عبارت‌دیگر محدب و مقعر بودن دامنه را 
نشان می‌دهد. انحنای پروفیل معرف اندازه تغییر شیب منحنی میزان» در طول مسیر جریان است و بنابراین انحنای 
پروفیل نشان دهنده شدت‌جریان آب و فرآیندهای حمل و رسوب‌گذاری می‌شود. به‌طوری‌که مقدار منفی این انحناء 
سطوح محدب (کوژ) و مقدار مثبت Of‏ سطوح مقعر GIS)‏ را نشان می‌دهد. دامنه‌های محدب فرسایش بیشتری 
نسبت به دامنه مقعر و هموار دارد (شکل AY‏ 

فاصله از آبراهه: منطقه‌ای که در مجاورت آبراهه‌ها قرار دارد به‌مراتب پتانسیل فرسایش بیشتری نسبت به مناطق 
دورتر خواهد داشت. لذا هرچه فاصله از آبراهه‌ها کمتر باشد. پتانسیل فرسایش بیشتر خواهد بود. در فاصله ۲۰۰- 
۰ متر بیشترین فرسایش رخ داده است (شکل AY‏ 

تراکم گسل: میزان تراکم گسل‌ها تأثیر بسزایی بر میزان خردشدگی و فرسایش سنگ‌ها خواهد داشست. بیشترین 
میزان حطر فرسایش خاک در تراکم گسل با مقدار تراکم از ۳۷۲۱ متر تا 1۲۷۸ متری تراکم گسل‌هادیده می‌شود 
(شکل ۲). 

فاصله از گسل: گسل‌ها باعث شکسته و خرد شدن سنگ‌ها می‌شوند. میزان فرسایش خاک در مناطق نزدیک به 
گسل‌ها بیشتر به وقوع می‌پیوندد. بیشترین میزان فرسایش خاک در فاصله کمتر از ۱۵۰۰ متری رخ داده است (شکل 
(Y‏ 


سال دهم بررسی پتانسیل مناطق مستعد فرسایش با مدل‌های .... yey‏ 


شاخص رطوبت توپوگرافی: شاخص رطوبت توپوگرافی به‌عنوان یک عامل مهم در ارزیابی پتانسیل فرسایش 
شناخته شده است. قدرت فرسایش رواناب ازنظر سرعت > OL‏ پتانسیل تخلیه و ظرفیت حمل رسوب از طریق 


شاخص رطوبت توپوگرافی به دست می‌آید. بیشترین میزان فرسایش در ۸۳۹/۲ به وقوع پیوسته است (شکل AY‏ 
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شکل ۲- بررسی عوامل مؤثر در وقوع فرسایش خاک 


3720000 


۲-۳- متغیرهای ورودی به مدل‌ها 


بدین معنا که به هر یک از مدل‌ها بر اساس معیارهای Sd‏ شده رتبه‌ای داده می‌شود و در پایان پس از مقایسه مجموع 


۱۰۹ جغرافیا و مخاطرات محیطی شمارة چهلم 


رتبه‌های به دست آمده مناسب‌ترین مدل برای تخمین بهترین معادله نفوذ انتخاب می‌شود. شیوه رتبه‌بندی این مدل‌ها 
بدین صورت است که در هر روش مدلی که کمترین میزان خطای MAE‏ را دارد. پایین‌ترین رتبه» یعنی رتبه و به 
مدل‌ها یا ترکیبات ورودی دیگر در هر یک از دسته‌های ورودی بر حسب میزان معیار خطاهای مذکور رتبه‌های ۲ 
الى آخر تعلق می‌گیرد؛ اما در مورد معیارهای ضریب همبستگی R)‏ مدل یا ترکیبی که بیشترین مقدار این معیارها را 
دارا است» رتبه ۱ احتصاص می‌یابد. همچنین به مدل‌هایی که میزان معیار خطای آن‌ها یکسان ALL‏ رتبه یکسان 
تعلق می گیرد. درنهایت پس از انتخاب مدل‌ها یا ترکیبات مناسب هر روش و مقایسه بین آن‌هاء روش برتر در منطقه 
موردطالعه انتخاب می‌گردد. بررسی مقایسه مدل‌های ICONA‏ ماشین بردار پشتییان» چاید و جنگل تصادفی نشان 
می‌دهد که می‌توان یک الگوی مشخصی برای منطقه مطالعاتی را با دقت مناسب برای برآورد فرسایش معرفی نمود و 
باید برای هر مدل» الگوی برتر را معرفی کرد. درنهایت پس از انتخاب مدل‌ها یا ترکیبات مناسب هر روش و مقایسه 
بین dl‏ روش برتر انتخاب می‌گردد. جدول ۲ ضریب همبستگی و خطای سناریوهای مختلف در چهار مدل را 
نشان می‌دهد. همچنین در این تحقیق تحلیل حساسیت به‌منظور شناسایی و تعیین درجه اهمیت هر یک از پارامترها 
در سری داده‌های ورودی بر عملکرد دینامیکی سیستم صورت گرفته است. به‌منظور آنالیز حساسیت در هر سری 
داده ورودی یک پارامتر حذف و مدل با همان داده‌های تست اجرا و مقادیر ضریب همبستگی و خطا تعیین گردید. 
با حذف هر پارامتر هر چه خطای مدل افزایش بیشتری نشان دهد Sly‏ تأثیر بیشتر آن پارامتر بر فرایند شبیه‌سازی 
است. جهت انجام تحلیل حساسیت برای سری داده‌های ورودی از مدل‌های CHAIDRF SVM-linear‏ و 


2۵ به دلیل دقیق‌تر بودن نتایج OF‏ استفاده شده است. 


جدول ۲- ضریب همبستگی و خطای سناریوهای مختلف در مدل ماشین بردار پشتیبان خطی. درخت تصمیم نوع 
جاید. مدل ICONA‏ و حنگل تصادفی در ایستگاه 


SVM-Linear Chaid Random Forest ICONA‏ سناریو 
M R MAE R MAE R MAE R MAE‏ 
M1 ۰/۸۹۳۱ Y/\00 ۰۹۹ ۱/۸۳۲ ۰ A ۰/۸۷۹ 7۹44‏ 
M2 ۰/۹۰۱ 744 ۰/۹۳ \/2£0 AAY \/00\ YAT ۳/۳۲‏ 
M3 ۰/۹۹۵ ۱۳۸/۸۵ AT \/0¥\ “/AVV ۱۳۲ ۰/۸۳/۹۵ ۳/۸۳۳"‏ 
M4 ۰/۷۱ ۱710۸ ۰۵ ۱/۳۳ ۰/۹۹۸ ۱۸/۳۱ AY ۱/۸٩۱‏ 
M5 ۰/۹۹۹ 1/04 ۰۹۹ ۳/۳۳" ۰/۹۵ AAKA) ۱ 744‏ 
M6 say ۱1:۶ AtA \/ow ۰۷۰ ۳/۳۳ ۰/۷۱ Y/ALO‏ 
M7 ۰/۸۹۹ ۱۱ ۰/۹/6۹ ۱/۷/۸۹ ۰/۹۳۱ YANAY NN) 1/77101‏ 
M8 ۰/۸۹٦ \NYO ۰/8۸ Y/AY\ ۰/4 3/0۲1 ۰/M4 1/40۳‏ 
M9 “/AVV 1/۸0۰ ۱ ۱/۳۵۳ AYY ۳/۳۲ ۰۳ YAYA‏ 
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SVM-Linear Chaid Random Forest ICONA‏ سناریو 
M R MAE R MAE R MAE R MAE‏ 
M10 46٥ + /MAV Nv) vava /40 ۱/۲ /ALO ۳/۰۳‏ 
M11 ۰۹۸ ۰,۸5۹ ۳۳ V/AV\ ۳ ۱/۳۳۳ ۹ 1/0۰۲‏ 
M12 ۰/4 ۱ iS) ۱/07 0 “AYA JAMA ۱۸۷۳۰۵‏ 


در پژوهش حاضر ۱۳ سناریو با ترکیب پارامترهای مختلف به‌عنوان ورودی مدل ایجاد شدند. نتایج این تحقیق 
نشان می‌دهد که در مدل ماشین بردار پشتیبان حطی بهترین سناریو M7‏ با بالاترین ضریب همبستگی ۰/۸۹۹ و 
کمترین مقدار خطا۴=۰/۷۱۱ M۸‏ مدل ICONA‏ بهترین سناریو ۷19 با بالاترین ضریب همبستگی ۰/۸۶۵ و کمترین 
مقدار خحطا MAE= VY)‏ مدل جنگل تصادفی بهترین سناریو 112 با بالاترین ضریب همبستگی ۰/۹۲۱ و کمترین 
مقدار خطا MAE=+/TWA‏ در مدل درخت تصمیم نوع چاید بهترین سناریو ۷10 با بالاترین ضریب همبستگی 
۷۱ و کمترین مقدار خطا ۷۸17-۱/۲۱۳ بوده است. درنهایت مدل جنگل تصادفی و ماشین بردار پشتیبان به‌عنوان 
بهترین الگوها با بیشترین ضریب همبستگی و کمترین خطا به‌عنوان الگوی برتر شناخته شدند. نتایج آنالیز حساسیت- 
ها نشان می‌دهد در مدل ماشین بردار پشتیبان حطی به ترتیب پارامترهای شیب کاربری اراضی, ارتفاع؛ تراکم شبکه 
آبراهه» خاک و در مدل جنگل تصادفی به ترتسی پارامترهای ارتفاع» شیب زمین شناسی. رطوبت توپوگرافی و 
کاربری اراضی به‌عنوان پارامترهای تأثیرگذار در مدل شناخته شدند. در مدل درخحت تصمیم نوع چاید پارامترهای 
جهت شیب. خاک ارتفاع و فاصله از آبراهه و در مدل ICONA‏ پارامترهای انحنا سطح» جهت شیب تراکم گسل و 
شیب به‌عنوان مهم‌ترین پارامترها؛ اهمیت نسبی بالایی در برآورد فرسایش خاک در منطقه موردمطالعه دارند. شکل 


(۳) پراکنش داده‌های مشاهداتی و برآورد شده توسط بهترین الگو از هر مدل نشان داده شده است. 


۸ جغرافیا و مخاطرات محیطی شمارة چهلم 


پراکنش مدل ICONA‏ پراکنش مدل تصادفی 
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شکل ۳- پراکنش داده‌های مشاهداتی هر مدل 


نتایج حاصل از ارزیابی فرسایش خاک بر اساس مدل‌های ICONA‏ ماشین بردار پشتیبان چاید و جنگل 

تصادفی در شکل‌های (۵) آورده شده است. نتایج به دست آمده از ارزیابی فرسایش خاک مشخص شد که در مدل 
ICONA‏ ۱۰ از سطح منطقه دارای فرسایش he‏ کم 7۱۳ فرسایش کم ۱۸/ فرسایش متوسط ۲۰/ فرسایش 
زیاد و ۸۳۹ فرسایش خیلی زیاد است. در ماشین بردار پشتیبان. ۱۲./ از سطح منطقه دارای فرسایش خیلی کم TNO‏ 
فرسایش کم ۸۱۳ فرسایش متوسط ۸۲۱ فرسایش زیاد و ۸۳۹ فرسایش خیلی زياد است. در مدل چاید ۸۳۲ از 
سطح منطقه دارای فرسایش خیلی کم» ۸/ فرسایش کم AYY‏ فرسایش متوسط. ۱۵ فرسایش زیاد و ۸۲۳ فرسایش 
خیلی زیاد است. در مدل جنگل تصادفی ۱1 از سطح منطقه دارای فرسایش خیلی کم 7۱۳ فرسایش کې ING‏ 
فرسایش متوسط. AYN‏ فرسایش زیاد و ۳۱/ فرسایش خیلی زیاد است (شکل AE‏ 
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شکل ۵- پهنه‌بندی فرسایش خاک بر اساس مدل‌های ICONA‏ ماشین بردار پشتیبان. جاید و حنگل تصادفی 


۶- جمع‌بندی 

در پژوهش حاضر برای ارزیابی be‏ فرسایش خاک در حوضه آبخیز گناباد از ۱۲ شاخحص Spe‏ در فرسایش 
خاک شامل شیب. سازند زمین شناسی» کاربری اراضی» خاک ارتفاع جهت شیب انحنا سطح» تراکم شبکه آبراهه 
فاصله از آبراهه» تراکم گسل» فاصله از گسل و شاحص رطوبت توپوگرافی Twi)‏ استفاده گردید. پارامترهای شیب 
سازند زمین‌شناسی» کاربری اراضیء خاک» فاصله از آبراهه و رطوبت توپوگرافی (Twi)‏ حساسیت بیشتری نسبت به 
فرسایش داشته و عوامل مؤثر بر فرسایش در این مناطق فعالیت بیشتری دارند. در منطقه موردمطالعه, به‌عصوص در 
اراضی فاقد پوشش گیاهی و اراضی مسکونی؛ فرسایش خاک نگران‌کننده است. اراضی آبرفتی جوان از نظر خاک 


۱۹۰ جغرافیا و مخاطرات محیطی شمارة چهلم 


زمین های زراعی استفاده می‌شود. این اراضی BE‏ محدودیت یا دارای محدودیت کم برای استفاده زراعت بوده و بايد 
با انجام عملیات صحیح زراعی از تخریب اراضی جلوگیری کرد. برای رسیدن به این منظور آموزش کشاورزان و 
برپایی دوره‌های ترویجی ضروری به نظر می‌رسد. خطر دیگری که این اراضی را تهدید می‌کند وجود فرسایش است 
که با انجام عملیات صحیح خاک‌ورزی در راستای کنترل فرسایش بایستی اقدام نمود. در ارضی کوهستانی مهم‌ترین 
مسئله وجود رخنمون سنگی و فرسایش آبی ناشی از تخریب پوشش گیاهی و یا استفاده نادرست از اراضی (زراعت 
دیم) و پستی‌وبلندی است. با توجه به حساسیت سازندهای حوضه به فرسایش بهترین راه جهت کنترل فرسایش 
احیاء و تقویت پوشش ALS‏ و تبدیل اراضی زراعی ديم به مرتع و یا کشت‌های جایگزین (زراعت‌هایی که عملیات 
خاک‌ورزی کمتری لازم دارند) است. با توجه به اينکه بخش عمده‌ای از زراعت‌های دیم در تیپ فلات‌ها و تراس- 
cl‏ فوقانی انجام می‌شود و با عنایت به اين‌که این رسوبات دارای حساسیت زیادی به فرسایش هستند و بایستی در 
استفاده از این اراضی مسائل حفاظت خاک را در نظر گرفت. در این اراضی زراعت های دیم بایستی به مرتع و با 
کشت‌های پایدار تبدیل شده و در مناطق شیبدار که پتانسیل فرسایش بیشتری وجود دارد باید تدابیر لازم اندیشیده 
شود. مناطق با شیب ۵/ حداکثر فرسایش‌پذیری مشاهده شده است که تحت کشت دیم قرار داشتند و اراضی فقیر 
است. افزایش شیب موجب می‌شود سرعت و حرکت هرز آب سطحی زیاد و مدت‌زمان تماس آب با خاک و 
درنتیجه نفوذ آب در خاک کم شود و همچنین موجب می‌شود که خاک بیشتر و سریع‌تر فرسایش یابد. ميزان تراکم 
گسل‌ها در قسمت‌های شرقی و جنوبی واقع شده‌اند و با افزایش فاصله از گسل فرسایش کمتری به وقوع می‌پیوندد. 
مقادیر فرسایش از قسمت‌های شرقی تا نیمه جنوبی روند کاهشی داشته و در قسمت‌های شمالی که ارتفاعات و 
بارش بیشتر است» فرسایش افزایش adh‏ است. ارتفاعات شمال شرق و شمال غرب به دلیل شیب زیاد بیشترین 
مقدار فرسایش را دارا می‌باشند. بررسی شبکه آب‌ها در منطقه نشان می‌دهد که شعبه‌های فرعی آب‌ها در مناطق 
شمالی پرآب‌تر و دائمی‌تر از شعبه‌های فرعی هستند. این در حالی است که تراکم شبکه آب‌ها در قسمت‌های 
غربی و جنوبی حوضه بیشتر از قسمت‌های دیگر آن است که به‌احتمال ناشی از نفوذناپذیری سنگ‌بستر زمین‌های 
این مناطق است. در مورد فاصله از رودخانه نیز با افزایش فاصله از رودخانه میزان پتانسیل فرسایش سطحی کاهش 
می‌یابد هرچقدر آبراهه‌ها به هم فشرده‌تر باشند. فرسایش هم زیادتر بوده است. ازآنجاکه بیشتر خاک‌های منطقه آنتی 
سول است و بافت OF‏ منفصل است و خاک فرسایش‌پذیری بیشتری دارد. شرایط و ویژگی‌های خاک و پوشش 
گیاهی بیشترین ارتباط را با عامل جهت شیب دارد که نهایتاً بر میزان پتانسیل فرسایش خاک در زمین دارد. بر اساس 
نتایج بیشترین مقدار فرسایش خاک در قسمت‌های مالی و جنوب غربی رخ داده است. از طریق خروجی 
مدل‌های ICONA‏ ماشین بردار پشتیبان tle‏ و جنگل تصادفی نقشه‌های حساسیت منطقه به فرسایش ang SLE‏ 
گردید. ارزیابی مدل‌ها با استفاده از ضریب همبستگی R)‏ و مقادیر میانگین قدرمطلق خطا (MAE)‏ صورت گرفت. 
بررسی نتایج مدل‌ها نشان داد که مدل‌های ماشین بردار پشتیبان 100۸4و جنگل تصادفی با سناریوی M7‏ 1۷9و 


سال دهم بررسی یتانسیل مناطق مستعد فرسایش با مدل‌های .... ۱۱۱ 


2 با بیشترین میزان همبستگی و کمترین خطای میانگین نسبت به مدل چاید به دلیل اینکه اساس OF‏ بر پایه قوانین 

منطقی بیان شده است از دقت بالایی برخوردار است. مقادیر مثبت سناریوها به معنای عملکرد ضعیف دیگر مدل‌ها 

نیست بلکه نمودی است که مدل‌سازی فرسایش تحت تأثیر آن است. ارزیابی و پیش‌بینی فرسایش خاک اطلاعات 

مفیدی را برای حفظ محیط‌زیست و حفاظت خاک ارائه می‌نماید. برای جلوگیری از تخریب زمین و آلودگی منابع 

آبی سطحی و زیرزمینی» پيشنهاد می‌شود که عملیات آبخیزداری و حفاظت خاک اعمال شود. نقشه‌های حاصل از 

پهنه‌بندی خطر فرسایش خاک می‌تواند مورداستفاده مهندسین منابع طبیعی و آبخیزداری قرار گیرد؛ زیرا با شناسایی 

پهنه‌های در معرض فرسایش خاک می‌توان با اقدامات مناسب از خساراتی که این پدیده بر محیط‌زیست. آلودگی 

منابع آبی سطحی و زیرزمینی و زمین‌های کشاورزی می‌گذارد. جلوگیری به عمل آورد و آن را کنترل کرد. 
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